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Abstract —  This  study  is  concerned  with  the  development 
of  an  integrated  multiple  compartment  model  of  the  iso¬ 
lated  cerebral  artery  in  rat.  The  smooth  muscle/arterial 
wall  complex  is  an  important  component  of  the  circulatory 
model  and  serves  as  an  “vasomotor  organ”,  which  provides 
the  myogenic  mechanism.  We  have  focused  on  this  myo¬ 
genic  mechanism  and  have  developed  a  model  of  the  elec- 
trophysiological,  contractile  and  mechanical  characteristics 
of  the  single  smooth  muscle  cell.  This  cell  model  is  used  to 
interrelate  the  topics  of  arterial  wall  stress,  changes  in  trans¬ 
membrane  potential,  intracellular  [Ca2+]i  concentration  and 
smooth  muscle  contraction.  Moreover,  the  small  cell  model 
is  embedded  in  a  larger  arterial  wall  model,  which  converts 
contractile  activity  into  changes  in  lumen  diameter.  The 
complete  model  is  used  to  provide  biophysically  based  ex¬ 
planations  of  the  myogenic  mechanisms  underlying  the  au¬ 
toregulation  of  cerebral  blood  flow. 
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I.  Introduction 

A  luminal  pressure  increase  in  cerebral  arterioles  and 
small  arteries  produces  an  initial  passive  distension  fol¬ 
lowed  by  a  sustained  constriction  due  to  active  force  de¬ 
velopment  in  vascular  smooth  muscle  cells.  This  process  is 
termed  the  myogenic  response  and  it  plays  a  critical  role 
in  maintaining  constant  blood  flow  and  pressure  for  micro- 
circulation. 

Several  models  have  been  developed  to  investigate  var¬ 
ious  aspects  of  vascular  smooth  muscle  function  includ¬ 
ing  Ca2+  dynamics,  contractile  kinetics  and  force  gener¬ 
ation  [1],  [2],  [3].  Models  have  also  been  developed  for 
isolated  cerebral  arteries  [4],  [5].  These  topics  are  usu¬ 
ally  studied  individually.  An  overall  modeling  approach, 
which  integrates  most  functional  compartments  in  different 
physiological  levels  including  membrane  electrophysiology 
of  smooth  muscle  cell,  smooth  muscle  contractile  mecha¬ 
nism  and  arterial  mechanical  properties,  is  needed  to  give 
a  complete  picture  and  better  understanding  of  this  physi¬ 
ological  system  behavior  with  regard  to  myogenic  response. 

In  the  current  study,  we  developed  a  mathematical  model 
to  describe  the  electrophysiological,  biochemical  and  me¬ 
chanical  processes  of  single  vascular  smooth  muscle  cell, 
and  hence  the  cellular  aspects  involved  in  the  mechanism 
of  myogenic  response.  The  single  cell  model  is  considered 
as  a  basic  functional  unit,  and  is  subsequently  coupled  to 
a  mechanical  model  of  arterial  vessel  which  describes  the 
elastance  and  pressure-volume  relationship  of  the  isolated 
cerebral  artery.  This  integrated  model  serves  as  a  gen¬ 
eral  framework  that  can  be  used  to  examine  physiological 
hypotheses,  and  also  provide  mechanistically  based  expla¬ 


nations  of  experimental  phenomena. 

II.  Model  Development 
A.  Model  of  Single  Smooth  Muscle  Cell 
A.l  Model  of  Cell  Electrophysiology 

Fig.  1  shows  the  diagram  of  the  model  for  smooth  muscle 
cell  electrophsyiology. 


A  Extracellular  space  . 


Fig.  1.  Cell  Electrophysiological  model  .  A.  The  model  describ¬ 
ing  membrane  currents  and  transmembrane  potential.  B.  The  fluid 
compartment  describing  ionic  dynamics  and  Ca2+  buffering,  includ¬ 
ing  Ca2+  handling  by  sarcoplasmic  reticulum  (SR). 

Fig.lA  represents  a  Hodgkin-Huxley  type  electrical 
equivalent  circuit  of  the  smooth  cell  membrane,  which  con¬ 
sists  of  a  whole-cell  membrane  capacitance  ( Cm )  shunted 
by  a  variety  of  resistive  transmembrane  channels  includ¬ 
ing  the  L-type  Ca2+  channel  (Ica,L),  mechanically  sensi¬ 
tive  channel  (Im),  Ca2+  activated  I\+  channel  (Ik,Co), 
as  well  as,  the  inward  rectifier  (Tat),  delayed  rectifier  (Ik) 
and  background  ( Ib )  currents.  Na+ /K+  and  Ca2+  pump 
currents  ( Ino,k  and  leap) and  the  Na+/Ca2+  exchanger 
current  ( iNaCa )  are  also  included.  Kirchhoff’s  current  law 
applied  to  the  equivalent  circuit  of  Fig.lA  yields  the  follow¬ 
ing  differential  equation  describing  changes  in  transmem¬ 
brane  potential  Vm : 

dVm  1 

— —  =  —  {ICa,L+lK+lK,Ca+lKi+lM+lNaCa+lNaK+ICaP+lB) 

at  C  m 

(1) 
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To  account  for  changes  in  the  ionic  concentrations  of 
Na+ ,  K+  and  Ca2+  with  voltage  and  time,  a  fluid  com¬ 
partment  model  (see  Fig. IB  )is  developed,  that  is  based  on 
cytosolic  ion  balances  for  Na+,  Ca2+  and  K+ .  Concen¬ 
trations  for  these  ions  are  assumed  to  be  constant  in  the 
larger  extracellular  bathing  medium.  The  material  balance 
for  Ca2+  accounts  for  membrane  Ca2+  influx  from  L-type 
Ca2+  channel  (major  pathway  of  Ca2+  entry),  buffering  by 
the  cytosolic  protein  calmodulin  (CM),  and  other  buffer¬ 
ing  media,  as  well  as,  a  secondary  ryanodine  sensitive  Ca2+ 
induced  Ca2+  release  (CICR)  mechanism  of  sarcoplasmic 
reticulum  (SR)  and  Ca2+  extrusion  via  sarcolemmal  and 
SR  membrane  Ca2+  pumps. 

This  model  is  based  on  various  experimental  studies  re¬ 
ported  in  the  literature  and  is  able  to  predict  responses 
to  a  variety  of  biophysical  processes,  most  importantly  the 
regulation  of  [ Ca2+]i . 

A. 2  Ca2+-CM  Dependent  Contractile  Kinetics 

The  interaction  between  actin  thin  filaments  and  myosin 
thick  filaments  provides  the  molecular  basis  for  muscle 
contraction  which  is  regulated  by  intracellular  calcium- 
calmodulin  complex  (CaCM).  This  complex  binds  to  and 
activates  the  enzyme  myosin  light  chain  kinase  (MLCK). 
Activated  MLCK  catalyzes  the  phosphorylation  of  myosin 
light  chain,  which  results  in  an  increase  in  actin-activated 
ATPase  activity  of  myosin  and  initiates  cross-bridge  cycling 
and  force  generation  [6],  [7],  [8]. 

We  employ  the  four-state  kinetic  model  developed  by  Hai 
and  Murphy  [2],  [9]  to  describe  myosin  phosphorylation  and 
latch  bridge  formation.  Fig.  2  shows  this  model,  which 
consists  of  four  fractional  species:  free  cross-bridges  (M), 
phosphorylated  cross-bridges  (Mp),  attached  phosphory- 
lated  cross-bridges  (AMp)  and  attached  dephosphorylated 
cross-bridges  i.e.,  latch  bridges  (AM).  The  latch  state  has 
been  found  to  be  unique  to  smooth  muscle,  and  it  endows 
the  smooth  muscle  cell  with  the  ability  to  maintain  force 
when  the  level  of  calcium-dependent  myosin  phosphoryla¬ 
tion  is  relatively  weak. 

The  population  of  phosphorylated  myosin  is  defined  as 
the  sum  Mp  +  AMp ,  whereas  attached  cross-bridges  are 
represented  by  the  sum  AMp  +  AM.  The  CaCM  depen¬ 
dence  of  myosin  phosphorylation  is  represented  through 
rate  constants  K\  and  Kq.  We  assume  Afl  =  Kq  which 
means  that  the  phosphorylation  process  has  same  rate  in 
changing  from  M  to  Mp,  as  from  AM  to  AMp.  Thus, 

I<i  =  K6  =  f(CaCM)  (2) 


A. 3  Smooth  Muscle  Cell  Mechanics 

To  study  the  active  mechanical  response  and  viscoelas¬ 
ticity  of  the  single  smooth  muscle  cell,  a  model  of  cell  me¬ 
chanics  is  developed  to  represent  the  coupling  between  the 
active  force  generation  of  myosin  kinetics  and  the  mechan¬ 
ical  property  of  the  myofilaments  and  the  cell  (Fig. 3). 
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Fig.  2.  Multi— state  kinetic  model  of  calcium  dependent 

myosin  phosphorylation  and  cross— bridge  formation 


The  elastic  response  originating  within  the  myofilaments 
of  the  cell  is  attributed  to  cross-bridge  stiffness  kx  in  se¬ 
ries  with  an  series  elastic  element  ks  having  an  exponential 
length-force  (L:F)  relationship  [10],  [11].  The  passive  elas¬ 
ticity  of  the  cell  is  modeled  as  a  parallel  element  kp ,  which  is 
characterized  by  nonlinear  exponential  (L:F)  relationship. 

The  characteristics  of  the  cross-bridges  is  described  ac¬ 
cording  to  the  force: 


_g(  la  lopt\ 2 

Fx  =  (kxi  ■  AMp  +  kx2  ■  AM)  ■  lx  ■  e  A  > 


(3) 


where  lx  is  the  extension  of  cross-bridge  and  lopt  is  the  op¬ 
timal  length  of  the  overlap  between  myosin  filament  and 
actin.  Both  latch  bridges  and  phosphorylated  attached 
cross  bridges  contribute  to  the  elasticity  of  this  compo¬ 
nent.  The  spring  constants  kx\  and  kx 2  in  equation  (3)  de¬ 
note  the  maximal  stiffness  which  can  be  achieved  by  latch 
bridges  and  phosphorylated  attached  cross-  bridges,  respec¬ 
tively.  AMp  and  AM  are  considered  as  the  outputs  of  the 
kinetic  contractile  model  described  in  previous  section,  and 
indicate  the  distribution  of  attached  cross-bridges  and  the 
state  of  interaction  between  myosin  and  actin  filaments. 

A  model  of  active  force  generation  is  developed,  based 
on  the  cross-bridge  kinetics  and  the  interaction  between 


Fig.  3.  Mechanical  model  for  smooth  muscle  cell  The  two 
branches  of  the  network  represent  the  mechanical  properties  of  the 
smooth  muscle  cell:  kp  denotes  the  component  for  passive  parallel 
elasticity,  whereas  kx  and  ks  denote  the  elasticity  of  attached  cross 
bridge  and  series  elastic  component,  respectively.  The  active  con¬ 
tractile  (length  la)  element  expresses  the  interaction  between  cross¬ 
bridges  and  actin  filaments  by  cross-bridge  cycling  and  sliding. 


the  cross-bridges  and  actin  filament.  The  force  developed 
by  this  element  is  described  as: 

Fa  =  [fAMp  ■  AMp  ■  (vx  +  la)  +  f  am  •  AM  ■  la\  ■  e  4  ‘°p*  ) 

(4) 

where  / AMp  and  f am  are  friction  constants  for  latch 
bridges  and  phosphorylatecl  cross  bridges,  respectively. 
The  assumptions  used  by  this  formula  are:  (a)  only  at¬ 
tached  phosphorylated  cross  bridges  have  cycling  behavior 
moving  actin  filament  with  a  linear  velocity  vx.  Note  that 
this  is  the  only  process  which  contributes  to  the  active  force 
generation.  The  latch  bridges  lose  cycling  capability  and 
remain  attached  to  actin  filament  changing  the  stiffness  of 
the  cell  (described  in  equation  (3));  and  (b)  both  species  of 
attached  cross  bridges,  phosphorylated  and  latch  bridges, 
can  move  along  the  actin  filament  with  a  velocity  la  in 
response  to  cell  stretch  and  release. 

B.  Vascular  Wall  Model 

We  also  consider  a  macroscopic  lumped  model  of  the  iso¬ 
lated  perfused  cerebral  artery,  in  terms  of  its  elastic  proper¬ 
ties  and  the  properties  of  the  perfusion  system  to  which  the 
vessel  is  attached  via  cannulation.  The  dynamic  response 
of  the  whole  vessel  to  a  pressure  pulse,  causes  changes  in 
vessel  diameter  that  are  mediated  by  the  stretch-activated 
response  of  individual  smooth  muscle  cells.  We  define  a 
lumped  elastance  function  E{t)  for  the  vessel  that  charac¬ 
terizes  the  pressure  -volume  relationship  for  this  distensible 
vessel.  The  smooth  muscle  cell  model  is  incorporated  into 
this  larger  vessel  model  and  provides  the  dynamic  changes 
involved  in  the  change  in  vessel  elastance  with  time.  Specif¬ 
ically,  we  develop  a  vascular  wall  model  based  on  a  cir¬ 
cumferential  distribution  of  smooth  muscle  cells  within  the 
arterial  wall.  This  is  formulated  as  a  model  of  vessel  elas¬ 
tance  E{t),  the  instantaneous  pressure- volume  relationship 
of  the  isolated  vessel,  which  is  modulated  by  the  activity 
of  smooth  muscle  in  the  wall.  Vessel  stiffness  and  the  hy¬ 
draulic  resistance  Ri  offered  to  blood  flow  are  altered  with 
changes  in  vessel  elastance  E(t).  This  leads  to  the  myo¬ 
genic  control  of  blood  flow  in  cerebral  vessel.  The  model 
diagram  is  shown  in  Fig. 4 

III.  Model  Simulations  and  Capabilities 

This  integrated  model  is  used  to  simulate  various  exper¬ 
imental  studies  on  both  the  single  isolated  smooth  muscle 
cell  and  the  isolated  arterial  vessel.  When  applied  to  study 
the  myogenic  mechanism  of  cerebral  artery,  the  model  is 
capable  of  mimicking  the  experimental  work  of  Knot  et 
al.  [12] .  Specifically  the: 

1.  Steady-state  relationship  between  intravascular  pressure 
and  membrane  potential; 

2.  Steady-state  relationship  between  intravascular  pressure 
and  arterial  wall  Ca2+  concentration; 

3.  Steady-state  relationship  between  intravascular  pressure 
and  arterial  wall  Ca2+  diameter; 


Fig.  4.  Diagram  for  the  vascular  model  of  isolated  cerebral 
vessel  and  testing  apparatus. 

4.  Membrane  potential  manipulation  by  [i^+]0  in  the  pres¬ 
surized  cerebral  artery; 

5.  Ca2+  dynamics  and  active  constriction  of  cerebral  ar¬ 
teries  in  response  to  an  elevated  intravascular  pressure. 

IV.  Summary 

This  integrated  model  provides  insights  summary  into 
the  cellular  mechanisms  involved  in  the  myogenic  autoreg¬ 
ulation  of  cerebral  blood  flow. 
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